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Abstract 

Taken by itself, the interference with the 7/ appears to strongly affect 
the amplitude of the transition 77 — ^ 37r. I point out that this effect is 
fictitious and also occurs in the mass spectrum of the pseudoscalars. 
Chiral symmetry implies that the same combination of effective cou- 
pling constants which determines the small deviations from the Gell- 
Mann-Okubo formula also specifies the symmetry breaking effects in 
the decay amplitude and thus ensures that these are small. 
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The decay 77 — > 37r is of particular interest, because it violates isospin 
symmetry. The electromagnetic interaction is known to produce only very 
small corrections [|TJ. Disregarding these, the transition amplitude is propor- 
tional to md — m u and thus represents a sensitive probe of the symmetry 
breaking generated by the quark masses. 

Since the quark mass term qmq = m u uu + ra^dd + m s ss breaks SU(3), 
it generates transitions between the octet and the singlet of pseudoscalar 
mesons. The consequences for the transition amplitude are discussed in the 
literature, but the results are contradictory: While a direct evaluation of the 
mixing effects [0 leads to the conclusion that the current algebra prediction 
is modified drastically, the chiral perturbation theory calculation to one loop 
[f|] yields the opposite result. The purpose of the present paper is to resolve 
this paradox. 

The direct calculation is based on the effective Lagrangian which describes 
the low energy structure of QCD in terms of a simultaneous expansion in 
powers of 1/N C , powers of the momenta p and powers of the quark mass 
matrix m. For a discussion of this framework to first nonleading order and 
references to the literature, see ref.||. The relevant effective field U(x) is an 
element of U(3) and includes the degrees of freedom of both the pseudoscalar 
octet and the singlet. The latter is described by the phase of the determinant, 
detU = e 1 ^ . Counting the three expansion parameters as small quantities of 
order l/N c = 0(5) , p = 0(y/5) , m = 0(5), the expansion starts with a 
contribution of order one, given by 

C eff = \F 2 (dJJ ] d»U) + \F 2 B(mU ] + mU) - + 0(5) , (1) 

where matrix traces are abbreviated with the symbol (...). The expres- 
sion involves three effective coupling constants: the pion decay constant 
F = 0(y/~Nc), the constant B = 0(1) which determines the magnitude of 
the quark condensate and the topological susceptibility r = 0(l). The rel- 
ative magnitude of the three leading contributions depends on the relative 
magnitude of the expansion parameters: the first involves two powers of 
momentum, the second is proportional to the quark mass matrix and the 
coefficient r of the third is smaller than F 2 or F 2 B by one power of 1/N C . 

Setting U = expi(p/F, the singlet field is given by the trace 4> = (ip)/F, 
so that the terms quadratic in ip are ^(d^ipd^ip) — ~B(m(p 2 ) — \r{ip) 2 /F 2 . 
For those fields which carry electric charge or strangeness, this expression is 
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diagonal and leads to the standard current algebra mass formulae, 

Ml + = {m u + m d )B , M 2 K+ = {m u + m s )B , M 2 K0 = {m d + m s )B . 

The states 7i°,f] and rj' undergo mixing. The mixing matrix is an element of 
0(3) and may thus be represented in terms of three angles, 9 m i, e, e'. The first 
arises from the mass difference between the strange and nonstrange quarks 
and breaks SU(3), while the other two are isospin breaking effects, driven by 
n^d — tti u . To first order in isospin breaking, the relation between the neutral 
components of the field, (p = ^3X3 + y?8^8 + V 9 9\/| an d the mass eigenstates 
7T°, 77, J]' is of the form 

<p 3 = 7r°-€rj-e'rj' (2) 
y? 8 = cos Q m > (i] + e tt°) + sin 9 m > (rj' + e' 7r°) 
V9 9 = - sin Q m , {i] + e n°) + cos Q m > {yj + e 7r°) . 

In the eigenvalues, the isospin breaking effects are of order (m^ — m u ) 2 . 
Neglecting these, the ir° is degenerate with tt ± , M 2 = (m u + m d )B. The 
remaining two eigenvalues involve a new scale, set by the topological suscepti- 
bility. Eliminating it, the diagonalization leads to two independent relations 
among the three quantities M v , M^, 9 m i, e.g. 

sin2^ = -|V2^^ , (3) 
Ml = |(4M| - Ml) + §V2 tg9 m , {M 2 K - Ml) . (4) 
The isospin breaking angles e, e' are proportional to the quark mass ratio 



\/3 (m d - m u ) , 

m = * 

4 (m s — m) 



e o = —— — , m = ^{m u + m d ) 



The coefficients of proportionality may be expressed in terms of 9. 

cos Q rm i — sin 

9 YjTj 1 \^2 



rfq' 1 



e = e cos 6 m i — (5) 

cos 9 m i + sin Q m i jy2 

cos 9 W +sin 

e = -2 e sin U w — -= . 

cos v m i — sin Vjrf V 2 
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Note that, in the counting of powers introduced above, M%, M^, M 2 , M"i, are 
treated as small quantities of order 5. According to eq.@, the mixing angle 
9 n > v is given by a ratio thereof and thus represents a quantity of order one. 
In this sense, the above formulae are valid to all orders in 6^. Numeri- 
cally, the ratio M 2 ,/M 2 is about equal to 3, indicating that the breaking of 
U(3) fl xU(3) L generated by the anomaly is larger than the breaking due to 
m s , by roughly this factor. The topological susceptibility, which describes 
the effects of the anomaly in the framework of the effective Lagrangian and is 
of order iV c °, is more important than the terms of order mN c , which account 
for the symmetry breaking generated by the quark masses. 

Consider now the decay r\ — > 7r + 7r~7r°. To calculate the corresponding 
transition amplitude with the effective Lagrangian in eq.(|I|), the expansion 
in powers of the field ip is needed to order if 4 . The first term yields a con- 
tribution proportional to ([d^ip, ^[d^ip, </?]). Since the remainder does not 
involve derivatives, the decay amplitude A involves at most two powers of 
momentum. Lorentz invariance and crossing symmetry then imply that A 
is of the form a + bs, where s = (p n + + p n -) 2 is the square of the center of 
mass energy of the charged pion pair. Performing the change of basis (@), one 
finds that b is given by — e/F 2 , where e is one of the mixing angles introduced 
above. The amplitude may thus be written as 

A = -e— (s-s A ) . 

The result is of the same structure as the current algebra prediction |J, 

A = -e j- 2 (s-lM 2 ) . 

As a consequence of the interference with the rj', the quark mass ratio e 
is replaced by the mixing angle e given in eq.(H). There is a corresponding 
change also in the value of the constant term, — > s A , but this term is 
inessential, for the following reason. In the limit m u ,m,i—>-0, the amplitude 
contains two Adler zeros, one at p w + = 0, the other at p n - = 0. For the 
above expression to have this property, the constant sa must tend to zero 
if m u , rrid are turned off. Hence an explicit evaluation would yield a result 
for sa proportional to M 2 : The interference with the rj' merely generates an 
SU(3) correction in the value | of the coefficient. Since contributions of order 
M 2 amount to small corrections, I will drop these in the following. 
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Using the observed values of M%, M|-, Ml, — M 2 as an input, the relation 
@ yields 9 m i ~ —22°, in reasonable agreement with what is found phe- 
nomenologically ||. Inserting this number, formula (^) gives e ~ 2eo- So, 
the net result of the above calculation is that mixing with the rj' increases 
the current algebra prediction for the amplitude by a factor of 2. 

I now compare this finding with the one loop result of chiral perturbation 
theory || . This calculation accounts for all effects of first nonleading order, in 
particular also for those due to r\r( mixing. It is based on SU(3) R x SU(3) L and 
hence only involves the degrees of freedom of the pseudoscalar octet. In this 
framework, the 77' only manifests itself indirectly, through its contributions 
to the effective coupling constants, like all other states which remain massive 
in the chiral limit, e.g. the p. 

Normalizing the amplitude with the kaon mass difference (e.m. self en- 
ergies removed) and with the pion matrix element of the axial current,^] the 
result is of the form 

^- (M ^:5 + W , M) , ( 6) 

3v 3 Fn 

where M(s, t, u) is a lengthy expression, which contains contributions gen- 
erated by the final state interaction, as well as symmetry breaking terms 
involving the effective coupling constants L 5 ,L 7 ,L 8 . The rf hides in the cou- 
pling constant L 7 , which also occurs if the mass of the 77 is calculated within 
the same framework. The explicit expression for the function M(s,t,u) con- 
tains this constant through a correction term which is proportional to the 
deviation from the Gell-Mann-Okubo formula and is denoted by 

_^ 4M 2 K -3M 2 -M 2 

^GMO — 



M 2 — M 2 

Dropping all other terms and disregarding contributions of order M 2 , the 
one loop result reduces to 



M(s ' t '" ) = 5Fm|(i + I^omo) + . . . = 4(MI 9 ! M|) (i + A OMO ) + . . . 



^ince the expression accounts for the corrections of order m, one needs to distinguish 
between the constant F in the effective Lagrangian and the observed decay constants 
F-kiFk, which differ from F through contributions of order to. 
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where I have used the identity (M*-M%)(1 + §A GMO ) = §(A/£-M*). The 
current algebra mass formulae quoted above show that, at leading order of 
the chiral perturbation series, the ratio (M^o — M^ + ) / {M\ — M%) is given by 
(m d —m u )/(m s —m) =4e /v / 3- Since the corresponding first order corrections 
do not involve the coupling constant L7, they are irrelevant in the present 
context and the same also holds for the difference between F v and F. 

With these simplifications, eq.(f|) reduces to A = — (e /F 2 ) s (1 + A GMO ): 
Up to small corrections of order Ml, the contribution from the symme- 
try breaking terms amounts to an overall renormalization of the amplitude, 
eo ~ >eo (1 + A GMO ). The experimental value of the deviation from the Gell- 
Mann-Okubo formula, A GMO = 0.22, shows that the modification is of rea- 
sonable size, confirming the general rule of thumb, according to which first 
order SU(3) breaking effects are typically of order 25%. The second order 
contributions, which the one loop formula neglects, are expected to be of the 
order of the square of this. Clearly, the outcome of the calculation described 
earlier is in flat contradiction with chiral perturbation theory. 

To identify the origin of the disagreement, I return to the earlier calcula- 
tion and express the angular factor occurring in eq.@ in terms of the masses 
of the particles. Solving the relation (Q) for tg6* w / and inserting the result, 
the angular factor becomes 



This is remarkable, because it shows that the expression for the mixing angle 
e may equally well be written as 



In this form, the result of the direct calculation differs from the corresponding 
term in the one loop prediction of chiral perturbation theory only by a factor 
of cos 6 V > V — 0.93, which represents a correction of order (m s — m) 2 and is be- 
yond the accuracy of the one loop result. I conclude that the two calculations 
are consistent with one another. In particular, it is incorrect to amalgamate 
the two by multiplying the one loop formula with the enhancement factor 
occurring in eq.(Q). 

The above expression for the angular factor shows that the result of the 
direct calculation is subject to an uncertainty comparable to the effect itself: 



cos^-sin^y/2 _ 2{2M 2 K - M 2 - Ml) 
cos 



1 + A, 



GMO 



e = e {1 + A GMO }cos^ 



(7) 
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Depending on whether one uses eq.(|3|) or eq.(|j), i.e. takes the observed 
values of Ml, — or as an input, the calculation yields e ~ 2eo or 
e ~ 1.2 eo, respectively. The problem arises from the fact that the mass 
formula (|4]) is not in good agreement with observation. If the second term 
is omitted, the relation reduces to the Gell-Mann-Okubo formula, which 
predicts M v ~ 566 MeV, slightly larger than what is observed. The second 
term indeed lowers the result, but the shift is much too large: While the Gell- 
Mann-Okubo prediction for M% only differs from the experimental value by 
7 %, the repulsion generated by mixing now yields a number which is too low 
by about 20 %. As pointed out in ref. H], early determinations of the mixing 
angle failed for precisely this reason: These were based on the assumption 
that the observed deviation from the Gell-Mann-Okubo formula is exclusively 
due to 777/ mixing and thus underestimated the magnitude of 9 V > V by about 
a factor of two. 

The above discrepancies do not indicate that the expansion of the effec- 
tive Lagrangian in powers of 1/N c ,p and m fails. Deviations of this order 
of magnitude are to be expected within a framework which only considers 
the leading term of the expansion. The effective Lagrangian in eq.(|l]) also 
predicts that Fx is equal to F n , while, experimentally, the two quantities 
differ by the factor 1.22. There is no reason why in the case of the masses, 
the corrections generated by the higher order contributions of the expansion 
should be smaller. 

The observed mass pattern is perfectly consistent with the assumption 
that the terms neglected in eq.(^) are small, but they definitely are differ- 
ent from zero. The main point here is that the same terms necessarily also 
affect the amplitude of the transition 77 — > 37r. The result of the chiral per- 
turbation theory calculation amounts to a low energy theorem: To order p 4 , 
the slope of the decay amplitude involves the same combination of couplings 
which determines the deviation from the Gell-Mann-Okubo formula. While 
the effective Lagrangian in eq. (|T]) only accounts for the coupling constant 
L 7 , which is related to r/r/' mixing, the one loop result receives significant 
contributions also from L5 and L$. In the framework of pole models 
these couplings are dominated by the exchange of scalar particles. Indeed, 
the particle data table shows that the mass of the lightest scalars is com- 
parable with M v >. These particles do not undergo mixing with pseudoscalar 
one-particle states, but with the ground state as well as with two-particle 
states. The corresponding effect in the square of the pseudoscalar masses is 
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also of order m 2 and is of opposite sign. It is suppressed by a relatively large 
energy denominator, because the masses of the scalars are large compared 
to those of the pseudoscalar octet, but the energy denominator is essentially 
the same as the one which suppresses the shift generated by 777/ mixing, so 
that the effects are of the same order of magnitude. 

This explains why the direct calculation does not yield a decent estimate 
for e, unless the result is written in the form ([?]), where it differs from the 
chiral perturbation theory result only by a factor of cos 6^. As noted above, 
this factor represents a second order correction of typical size. I conclude 
that there is no indication for the symmetry breaking effects of higher order 
to be unusually large - the factor cos 9 V ' V may be taken as an estimate for 
the uncertainties in the decay amplitude due to these. 

The current algebra prediction for the decay 77 — > 37r also receives cor- 
rections from a quite different source: final state interactions. These are 
generated predominantly by two-particle branch cuts and are responsible 
for the bulk of the one loop corrections. The corresponding higher order 
contributions may be worked out from unitarity, using dispersion relations 
and relying on chiral perturbation theory only to determine the subtraction 
constants || [T^] . Viewed in this perspective, the above discussion implies 



that the one loop predictions for the subtraction constants are trustworthy: 
These account for all symmetry breaking effects of first nonleading order, in 
particular for those due to 7777' mixing, and there is no indication for large 
corrections from higher orders. 

Note that these statements need not hold for radiative transitions, which 
may well be distorted by the pole due to //'-exchange @, |TlJ. The 3n channel 
is special, because the transition amplitude is determined by the effective 
Lagrangian of the strong interaction — this is why it is firmly tied to the 
mass spectrum of the pseudoscalars. 
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